Mon. Not. R. Astron. Soc. 000. HHTol (20021 Printed 2 February 2008 (MN WF&L style file v2.2) 



The complex structure of low luminosity active galactic 
nuclei : NGC 4579 



Marcella Contini * 

School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 



Accepted . Received ; in original form 



ABSTRACT 

We have modelled the low luminosity AGN NGC 4579 by explaining both the 
continuum and the line spectra observed with different apertures. It was found that 
the nuclear emission is dominated by an AGN such that the flux from the active centre 
(AC) is relatively low compared with that of the NLR of Seyfert galaxies. However, 
the contribution of a young starburst cannot be neglected, as well as that of shock 
dominated clouds with velocities of 100, 300, and 500 kins" 1 . A small contribution 
from an older starburst with an age of 4.5 Myr, probably located in the external 
nuclear region is also found. HII regions appear in the the extended regions (~ 1 
kpc), where radiation and shock dominated clouds with V s =100 kms _1 prevail. The 
continuum SED of NGC 4579 is characterized by the strong flux from an old stellar 
population. Emissions in the radio range show synchrotron radiation from the base of 
the jet outflowing from the accretion disc within 0.1 pc from the active centre. Radio 
emission within intermediate distances (10-20 pc) is explained by the bremsstrahlung 
from gas downstream of low velocity shocks (V s =100 kms -1 ) reached by a relatively 
low radiation flux from the AC. In extended regions (> 100 pc) the radio emission 
is synchrotron radiation created by the Fermi mechanism at the shock front. The 
shocks are created by collision of clouds with the jet. All types of emissions observed 
at different radius from the centre can be reconciled with the presence of the jet. 

Key words: galaxies: active-galaxies: nuclei-galaxies individual: NGC4579. 



1 INTRODUCTION 

The analysis of the line and continuum spectra of active 
galactic nuclei (AGN) developed in the last years, revealed 
that radiation in the different spectral ranges depends on dif- 
ferent excitation mechanisms of gas and dust. It was found, 
in particular (e.g. Contini & Viegas 2000, Contini, Viegas 
& Prieto 2004) that both the photoionizing flux from the 
active centre (AC) and collisional processes in supersonic 
velocity fields contribute to ionization and heating of the 
emitting gas (and dust) in the narrow emission-line region 
(NLR) with different relative importance. 

In low luminosity objects, such as LLAGNs (low lu- 
minosity AGN) and LINERs (low ionization nuclear emis- 
sion regions) starbursts as well as an active nucleus can be 
present with comparable importance. Anyway, winds and 
jets are often observed, leading to a prominent role of shocks. 
The complex nature of LINERs was revealed by modelling 
the most significant line ratios in the UV-optical-IR spectra 
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(Contini 1997). Moreover, modelling the X-ray -IR contin- 
uum correlation in a large sample of AGN, Contini, Viegas 
& Campos (2003) found that LLAGNs are in the low lumi- 
nosity tail of AGN. This is explained by a relative low flux 
from the AC, roughly two orders of magnitude lower than 
that found in the NLR of Seyfert 2 and Seyfert 1 galaxies, 

In this paper we analyse the spectra of NGC 4579 which 
is classified as a LINER 1.9 galaxy with many AGN charac- 
teristics by Pogge et al. (2000) and as a Seyfert 1.9 by Ulves- 
tad & Ho (2001). The morphological type of NGC 4579 is 
SBb. The galaxy is located at a distance of 20 Mpc (cz=1502 
kms -1 ) (Gonzalez-Delgado & Perez, 1996). Hubble Space 
Telescope (HST) images at 2300 A (Maoz et al. 1995) show 
emission resolved into two components, a bright core and 
a weaker source with a flux corresponding to 15 % of that 
of the bright core. The weaker source is extended 0.1 X 0.2 
arcsec 2 and is located 0.2 arcsec from the bright core at 
PA=72°. 

From ground-based observations the galaxy was previ- 
ously known to contain a type 1 nucleus revealed by a broad 
Ha line (FWHM of ~ 2300 kms -1 ). New spectra obtained 
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by HST 0.2" slit reveal an Ha component with FWZI of 
~ 18,000 kms _1 (Barth et al, 2001) with "shoulders" on the 
red and blue sides. About 31 low-luminosity HII regions have 
been detected and measured in the disc, located along the 
dust lanes of the spiral arms, at distances between 3.3 and 
7.5 kpc from the nucleus (Gonzalez-Delgado & Perez). 

NGC 4579 line spectra have been observed with dif- 
ferent apertures in the different frequency domains by 
Gonzalez-Delgado & Perez (1996) within 1.5" in the nuclear 
region and within 20" in regions A and B, and previously by 
Keel (1983) within 4.7", and by Stauffer (1982) within 4.0". 

In this paper we address the line spectra observed by 
Gonzalez-Delgado & Perez with small and large apertures. 
The latter include some information about the surround- 
ings of the nucleus. The earlier data by Keel and Stauffer 
obtained with intermediate apertures will be considered only 
to check the consistency of modelling. In fact, the large aper- 
tures integrate on all the conditions found in the nucleus, the 
values from the small aperture and the surroundings. In the 
case of this galaxy they are rather complex and includes ev- 
ery thing from a nuclear stellar bar, a starforming ring, an 
the old galaxy population, etc. 

Indeed, the low [OIII] 5007/4363 observed line ratio 
(Sect. 2.2) indicates relatively high temperatures in the 
emitting gas, because very high densities (> 1000 cm~ 3 ) 
downstream are excluded by [SII] 6717/6730 ^ 1. The radi- 
ation flux from the AC cannot heat the gas of the NLR to 
temperatures > 2 10 4 K, which, on the other hand, are eas- 
ily reached by collisional heating downstream of the shock 
front. Moreover, the [Oil] 3727+/H/3 (the + symbol indi- 
cates that the doublet intensities are summed up) and [OI] 
6300+/H/3 line ratios are relatively high. So pure photoion- 
ization models (cfr. the code CLOUDY) are less adapted to 
model NGC 4579. The SUMA code (Viegas & Contini 1994, 
Contini & Viegas 2001a, b and references therein) which cal- 
culates the spectra emitted from gas and dust accounting for 
the coupled effects of shocks and of an external photoioniz- 
ing source, is adopted throughout this work. 

The spectral energy distribution (SED) of the contin- 
uum is calculated consistently with the line spectra to con- 
strain the models. Particularly, the SED of NGC 4579 con- 
tinuum is characterised by an unusually high contribution 
of the old stellar population, higher by a factor of ~ 10 than 
in Seyfert 2 and Seyfert 1 galaxies. 

The observations in the radio by Anderson, Ulvestad, & 
Ho (2003) provide the fluxes within 10 4 Schwarzchild radii 
of the black hole, i.e. within 0.1 pc from the nucleus for a 
sample of LLAGN (NGC 3147, NGC 4168, NGC 4203, NGC 
4235, and NGC 4450). NGC 4579 was actually selected from 
this sample because radio data are available also on larger 
scales. It was suggested by Ulvestad & Ho (2001) that jets 
are responsible of the subparsec radio emission. This may 
originate in synchrotron radiation from discrete plasma com- 
ponents or from the "base" of a continuous jet ejected from 
the central engine, which becomes optically thin on larger 
scales (Nagar, Wilson, & Falcke 2001). Anderson et al. ex- 
plain the observed spectral slope by compact jets that can 
accomodate the average spectral index, the relatively high ra- 
dio luminosity, and the unresolved appearance. This applies 
only if the jets are fairly close to our line of sight. Anderson 
et al. conclude that, as the galaxies show a broad line com- 
ponent, they have to be in the right direction. This suggests 



that synchrotron emission with different spectral indices can 
be observed at different distances from the AC, and that 
shock fronts are at work. The shocks are created by colli- 
sion of jets with the interstellar medium (ISM). The in-situ 
particle acceleration associated with shocks could lead then 
to extended regions of radio emission with suitable spectral 
indices. Moreover, if the jets were disrupted by interaction 
with the ISM, the emission line and continuum spectra from 
the NLR clouds should show the characteristic features of 
shocks, such as relatively high low ionization (e.g. [Oil], [Nil], 
and [SII]) and neutral (e.g. [OI]) line ratios to H/3, and high 
[OIII] 4363 /[OIII] 5007+. Moreover, soft X-ray emission 
and dust reradiation in the IR could be explained, as well 
as faint diffuse radio emission (Middleberg et al. 2004). The 
signature of shocks is well recognizable in the radio range of 
the continuum because the power-law type of synchrotron 
radiation has a different slope than that of bremsstrahlung 
from cool gas (Contini & Viegas 2000). 

In this paper, we will investigate whether radio emis- 
sions within different radii from the AC on large scales are 
consistent with the jet hypothesis by showing that the dif- 
ferent types of emissions within different regions can be ex- 
plained by the shocks. Actually, shocks in the NLR of AGN 
are investigated by modelling the line and continuum spec- 
tra. In order to unravel the different aspects of the NGC 
4579 complex (e.g. the AGN, the starburst, and the HII re- 
gion), we will model NGC 4579 first through the line spectra 
(Sect. 2). We will calculate the different components of the 
continuum SED and compare them with the data observed 
with different apertures in Sect. 3. Discussion and conclu- 
sions follow in Sect. 4. 



2 THE LINE SPECTRA 

The spectra in the nuclear region observed by Gonzalez- 
Delgado & Perez (1996) within 1.5" at the 4.2 m William 
Herschel Telescope in La Palma during 1988 March, and 
at the 1 m Jacobus Kapteyn Telescope during 1992 May 
show a number of line ratios large enough for modelling. The 
spectra observed in the circumnuclear region (< 20") contain 
a small number of lines, and are not reddening corrected. 
They can however provide some precious results. The line 
spectra observed by Keel (1983) at 4.7" and by Stauffer 
(1982) at 4.0" at the Lick 1 m Anna Nickel telescope are 
relatively poor in number of lines. 

Pogge et al (2000) in their Table 3 give the ratio of 
[OIII]/Hq in the nuclear (0.23", 18pc) region and in the 
circumnuclear (2.5", 194pc) regions, which correspond to 
[OIII]/Hq=0.37 and 0.06 , respectively. However, they can- 
not be used for the modelling because they refer to "band" 
fluxes, not converted to emission in a particular line by cor- 
recting for the filter transmission of other lines in the band- 
pass. 

2.1 Modelling procedure 

We consider clouds moving outwards from the centre. The 
inner edge of the clouds is illuminated by the flux from the 
photoionizing source (AGN, starburst, HII region, etc) while 
the outer edge defines the shock front. 

As a first guess, we compare the data with models from 
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the grids calculated by the SUMA code (Contini & Viegas 
2001a,b, thereafter CVOla and CVOlb). The models pre- 
sented in (CVOla) correspond to the AGN with power-law 
flux intensities ranging from zero in shock dominated (SD) 
clouds up to a maximum of log Fh~ 13 which was found 
in radiation dominated (RD) clouds by previous modelling 
of AGN (e.g. Contini & Viegas 2000). Models presented in 
CVOlb refer to the starbursts where the ionizing flux cor- 
responds to a stellar cluster. Models for HII regions where 
black body radiation reaches the inner edge of the cloud are 
also presented in CVOlb. All the models are composite, i.e. 
they account also for the shocks, as suggested by the rela- 
tively strong low ionization level lines (Contini 1997), by soft 
X-ray, as well as by infrared luminosities (Contini, Viegas, 
& Campos 2003). 

The input parameters of SUMA relative to the shock are 
the following : the shock velocity, V s , the pre-shock density, 
no and the pre-shock magnetic field, Bo. Those relative to 
the flux are the power law flux, Fh, in units of photons cm~ 2 
s" 1 eV _1 at 1 Ryd, the spectral indices in the UV and X-ray 
ranges. a uv = -1.5 and a x = -0.4 are adopted, respectively. 
In models representing starbursts, t is the starburst age, and 
for HII regions T» is the color temperature of the stars. U 
is the ionization parameter in both cases. For all types of 
models D is the geometrical thickness of the cloud, and d/g 
the dust-to-gas ratio. A preshock magnetic field B = 10~ 4 
gauss and cosmic relative abundances (Allen 1973) are used 
for all the models. 

The intensity of each line depends on the physical con- 
ditions of the emitting gas, namely, on the distribution of 
the temperature, density, and fractional abundance of the 
corresponding ion downstream of the shock front, and on 
the relative abundance of the elements. 

In our previous modelling of the NLR of AGN (e.g. 
Contini et al. 2003 and references therein) it was found 
that multi-cloud models explain the observed spectra. These 
models result from the weighted sum of models representing 
single clouds each characterized mainly by the shock ve- 
locity, the gas density, and the radiation flux intensity. For 
LINERs, in particular, it was found that the single-cloud 
models account for all types of ionizing and heating mecha- 
nisms (Contini 1997), e.g. the power-law radiation from the 
active nucleus, the radiation from starburts, the black body 
radiation from stars with a certain colour temperature, and 
shocks created by collision of clouds in the ISM and in super- 
nova remnants (SNR). Therefore multi-cloud models will be 
adopted in the present modelling. The relative importance 
of the different types which best fit the observed spectra will 
shed some light on the characteristics of the galaxy on scales 
larger than the nuclear ones. 

2.2 Single-cloud models : the physical conditions 
of the emitting gas 

Single-cloud models are shown in Table 1. Some models that 
do not contribute to the best fit are also included (e.g. M5 
and M10) in order to understand the dependence of the line 
ratios on the different input parameters. 

Interestingly, the models which roughly fit the line ra- 
tios belong to all of the different types which are considered 
in the interpretation of LLAGN spectra (Contini 1997, Con- 
tini et al. 2003) : SD models with a relatively low (V s =100 
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Figure 1. [OII]/H/3 versus [OI]/H/3 calculated by the grids. Black 
triangles represent models calculated by a power-law radiation 
flux; open triangles refer to shock dominated models; open squares 
represent models calculated by black body, and black squares refer 
to starburst. The vertical solid line indicates [OI]/H/3= 0.75. 



kms Hn Ml) and high (V s =500 kms Hn M3) shock ve- 
locities, RD composite models corresponding to a relatively 
low V s and a low radiation flux from the AC reaching the 
cloud (M4, M5, M6, M7, and M8), composite models rela- 
tive to a young starburst (M9 and M10) and to a starburst 
with an age of 4.5 Myr (Mil), and a composite model cor- 
responding to an HII region (M12). Both radiation-bound 
and matter-bound models appear in Table 1. 

Cid Fernandez et al (2004) and Gonzalez-Delgado et al. 
(2004) claim that the relative importance between collisional 
phenomena and the photoionizing radiation which originates 
from the AC and/or from starbursts, or even from HII re- 
gions, can be deduced from the line ratios. Very recently, 
analysing new ground observations of LLAGN and data by 
HST, respectively, they distinguished between objects with 
high [OI] 6300/Hq (>0.25) and low [OI] 6300/Ha relatively 
to the stellar population age and to the presence of an active 
nucleus. The data collected from the grids are given in Fig. 
1 . It can be seen that the data referring to starburst models 
show [OI]/H/3< 0.75 ([OI]/Ha< 0.25, adopting Ha/H/3=3), 
out of few models calculated by V s =100 kms -1 , U=0.01, an 
age of t= 3.3-4.5 Myr, and D=3 pc. In other words, these 
clouds are more likely located in the surroundings of the 
starburst, and they are merging with the ISM. All the other 
models showing [OI]/H/3> 0.75 refer to composite models 
with a power-law flux, shock dominated models, and mod- 
els for HII regions. 

The spectra calculated by SD models with V s =300 
kms~ 1 (M2) and V s =500 kms _1 (M3) and a preshock den- 
sity no =300 cm" , are chosen following our previous results 
that the velocity gradient increasing towards the AC is ac- 
companied by a density gradient, and because [NIJ/H/3 are 
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Table 1. Single-cloud models 
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3.0 
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.0 
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.0 
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11.0 


.6 


.0 


2.0 
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.0 


.0 


.0 
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.0 
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Regions A and B (MC2) : 
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.0 
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.0 


.0 


.0 


.0 
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% [OI] 6300+ 
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.0 


.0 
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.0 


.0 
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.0 


.0 


.0 


.0 
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% H/3 
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.0 


.0 


27.0 


.0 


.0 
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.0 


.0 


.0 
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1 SD=shock dominated; pl=power-law; sb=starburst; MI: HII region 

2 in 10 -3 erg cm -2 s -1 

3 in photons cm~ 2 s _1 cV _1 at 1 Ryd 



relatively high and [SIIIJ/H/3 relatively low. Model M4 ac- 
counts for a relatively low ionization flux from the AC that 
is appropriate to LINERs. The composite model which ac- 
counts for a starburst at early age and a small ionization 
parameter (M9) shows a relatively low V s . As was pointed 
out by Viegas, Contini & Contini (1999) and confirmed by 
Contini & Contini (2003) low shock velocities are character- 
istic of young starbursts, which had not yet the time to de- 
velop supernovae. The only model corresponding to an older 
starburst (Mil) which could be considered corresponds to 
a low V s . Most of the line ratios to H/3 are higher than 
observed, out of [OIII] 4363 and [NI] 5200+ . Finally, the 
spectrum emitted from gas illuminated by a HII region cor- 
responding to a temperature of 5 10 4 K and U=0.1 (M12) 
is characterised by low [Oil] /H/3 line ratios. 

The single-cloud models fit the data very roughly. Yet 



they can lead to a better fit when summed up adopting 
relative weights. 

2.3 Multi-cloud models : comparison with the 
data 

Modelling the NLR of different galaxies (Seyfert types 1 and 
2, NLS1, luminous infrared starburst galaxies, LINERs, etc) 
we have learned that the spectra are generally explained by 
multi-cloud models. Models MCI and MC2 which appear 
in Table 2 result from the weighted averages of single-cloud 
models which better fit the line ratios. The weights corre- 
sponding to the different single-cloud models are given in 
the rows 13 and 14 from the bottom of Table 1. The relative 
weights compensate for the H/3 absolute flux intensities and 
indicate roughly the covering factors of each type of clouds. 
So e.g. the relative weight of model M12 is low, but this does 
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Table 2. Comparison of calculated with observed line ratios 
(H/3=l) 
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1.8+ 0.2 
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2.25 


2.49 


1.3 


[SII] 6730 


1.8+ 0.2 


2.2 


1.89 


1.86 


1.1 


[Oil] 7325+ 


0.61 


0.8 








[SIII]6069+ 


0.29+ 0.05 


1.8 









1 Gonzalcz-Dclgado & Perez (1996), nuc. (1.5"), corrected with 
C(H/3)=0.97 

2 Gonzalez-Delgado & Perez (1996), A « 20"), not corrected 

3 Gonzalez-Delgado & Perez (1996), B « 20"), not corrected 



not mean that HII regions do not contribute somewhat to 
the NGC 4579 spectra. 

In Table 2 we compare the observed line ratios with 
model results. The results show that the N/H relative abun- 
dance is higher than cosmic (9.1 1CT 5 ) by a factor < 2 
within 1.5" and by a factor > 3 in the extended region. 
The [OIII] 4363/5007 line ratios are underpredicted by a 
factor of 3, while [NIJ/H/3 and [SIII]/H/3 are underpredicted 
and overpredicted by factors of ~ 6, respectively in the nu- 
clear region. The fit, in fact could be improved by models 
calculated purposely for NGC 4579 and/or by a larger num- 
ber of single-cloud models. The average spectra correspond- 
ing to the extended regions A and B (MC2) overpredict 
the [OI]/H/3 line ratio and underpredict the [SII]/H/3 ratios. 
Generally, a lower [OI] flux corresponds to a lower [SII] flux, 
because the first ionization potential of S is lower than that 
of H and, therefore, of O. So, we conclude that also the S/H 
relative abundance should be about twice the cosmic one 
(1.6 10~ 5 ) in the extended region. The number of the ob- 
served line ratios is small and does not constrain strongly the 
averaged model. Nevertheless, the present analysis, which is 
carried out only by models from the grids, shows that the 
choice of models presented in the grids is sensible. 

The results show that in NGC 4579 the AGN coexists 
with starbursts and HII regions, as is generally found for 
LLAGNs. The relative importance of the different types is 
determined by the relative weights adopted summing up the 
single-cloud models. 

In the bottom of Table 1 the contribution in percent of 
each model to each of some significant lines ([Oil] 3727+, 
[OIII] 4363, [OIII] 5007+ , [OI] 6360+ , and H/3) is given for 
the nuclear region and in regions A and B. 

It can be noticed that the largest contribution (~ 60 %) 
to the nuclear spectra comes from RD clouds corresponding 
to the highest flux (log Fh= 10) and relatively low velocities 
(M8). SD clouds corresponding to V s =300 kms~ 1 (M2) also 
contribute by factors > 8.6%, with maxima of ~ 80 % to 
the [OIII] 4363 line. A contribution < 5% to all the lines 



comes from SD clouds with V s =500 kms _1 (M3). Interest- 
ingly, a young starburst (M9) contributes by less than 13 % 
with maxima of ~ 12 % to [SII] 6717 and H/3. Moreover, a 
starburst with an age of 4.5 Myr (Mil) contributes to all 
the lines. 

Shock velocities of 100 kms _1 in the young starburst are 
reasonable, because the starburst had not yet the time to 
develop SN. However, the low V s (100 kms -1 ) and n (100 
cm~ 3 ) in the older starburst may indicate that the starburst 
is in the outskirts of nuclear region where the outward ve- 
locities are generally low, while the young starburst is most 
probably correlated with the jet within the nuclear region. 
A similar situation was found in the Seyfert 2 galaxy NGC 
7130 (Contini et al. 2002) where a young starburst was re- 
vealed by modelling the spectra in the nuclear region, even 
if the effect of the AGN dominates, while the circumnuclear 
regions show spectra corresponding to stars with lower tem- 
peratures. 

Regions A and B show a strong contribution from SD 
clouds (Ml) and RD clouds (M4, M7) reached by a rather 
low flux, and having V s =100 kms" . HII regions corre- 
sponding to a temperature of 5 10 4 K (M12) are also present. 

It is clear that the HII regions are located in the circum- 
nuclear regions of NGC 4579, while the AGN and the young 
starburst coexist in the nuclear region. The distances of re- 
gions A and B from the nucleus are within 1 kpc (Gonzales- 
Delgado & Perez 1996, Fig. 4a). The modelling shows that 
high V s do not appear in these regions. 

Indeed, the classification (Cid Fernandes et al, 
Gonzalez-Delgado et al 2004) between LINERs and Tran- 
sition Objects considering the [OI]/Hcv line ratio sets NGC 
4579 among LINERs and is consistent with the presence of 
an AGN nucleus and of an old stellar population. 

Finally, the spectra observed within 4.7" and 4.0" by 
Keel (1983) and Stauffer (1982), respectively, hint about the 
conditions beyond the nuclear region. The spectra show sig- 
nificant contributions from HII regions (M12), SD clouds 
with V s =500 kms _1 (M3), and of an AGN with log F h =9 
(M4). It is typical of the NGC 4579 continuum SED that the 
old star population strongly dominates in the optical range 
(see Sect. 3). These stars do not affect the spectra because 
their temperatures are too low to heat and ionize the sur- 
rounding gas. The signature of old stars in the spectra could 
be revealed by collisional events in supernova winds leading 
to high velocity shocks which contribute to the heating and 
ionization of the gas in the clouds. 



3 THE CONTINUUM SED 

The modelling of line and continuum spectra must be cross- 
checked until the best fitting model is found for both. 

The bremsstrahlung emitted from gas in the down- 
stream region of the NLR clouds, generally explains the ob- 
served SED in a wide frequency range (10 8 - 10 18 Hz). In the 
infrared range, however, dust emission dominates, while the 
old stellar population background is observed in the opti- 
cal - near-IR. Synchrotron radiation generated by the Fermi 
mechanism at the shock front is often observed in the radio 
range of AGN. 

In Figs. 2 we present the SED of the continuum of NGC 
4579 in the different spectral ranges, while Fig. 3 shows the 
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Figure 2. The SED of the continuum. Top-left panel: the bremsstrahlung corresponding to the different models. M4 (short-dashed 
lines), M3 (solid lines), M8 (short-dash-long-dashcd lines) model corresponding to V s =1000 km s _1 (dotted lines), model with V s =100 
kms -1 , no=100 cm~ 3 (long-dashed lines); top-right : dust emission (same notation as for top-left); bottom left: black-body radiation 
from the old population stars; bottom right : the slopes in the radio range: synchrotron radiation from the Fermi mechanism (solid line), 
bremsstrahlung and dust emission from models M4 and M8. 



SED of the multi-cloud spectrum in a large frequency range. 
Open triangles represent data from the NED/IPAC Extra- 
galactic Database) (Fabbiano et al. 1992, Kinney et al. 1993, 
Kodaira et al. 1990, De Voucouleurs et al. 1991, Zwicky et 
al. 1961, Spinoglio et al. 1995, Aaronson 1977, Jarrett et 
al. 2003, Rousscl et al. 2001, Scoville et al. 1983, Soifer et 
al. 1989, Moshir et al. 1990, Tuffs et al. 2002, Whiteoak 
1970, Becker et al. 1991, Gregory & Condon, 1991, Dressel 
& Condon 1978, Heeschen & Wade 1964), where the data 



in the soft X-ray range are upper-limit. Open squares show 
the data from observations in the near-mid-IR by Alonso- 
Herrero et al (2003). Black triangles represent the data of 
Anderson, Ulvestad, & Ho (2003) from 1.7 to 43 GHz us- 
ing the Very Long Baseline Array, while the black squares 
show the data by Hummel (1980) obtained at the West- 
erbork Synthesis Radio Telescope, and the data by Crane 
(1977). Notice that Hummel presents two data at 1.5 GHz, 
namely, the flux density of the central source (< 21", i.e. 
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< 1.6 kpc) and the total flux density ol the galaxy with a 
radio extent of 3.9 arcmin along the major axis. The fluxes 
from Crane 1977 at 2.7 and 8.1 GHz with size < 3" (about 
230 pc, adopting a distance to Earth of 16 Mpc as given 
by Roussel ct al 2001) are also given by Hummel (1980). 
The data from the NED by Heeschen & Wade (1969) broad 
band measurement, as well as the total flux at 1.5 GHz by 
Hummel account for radiation from a more extended region. 

3.1 Bremsstrahlung 

Consistent modelling of line (Sect. 2) and continuum spec- 
tra (Contini, Viegas & Prieto 2004) of NGC 4579 shows that 
the SED is explained by a multi-cloud model which accounts 
for SD clouds with shock velocities of 500 km s -1 (solid lines) 
and RD clouds with a low V s (100 kms -1 ) reached by a rela- 
tively low photoionization flux from the AC (log F h =9) rep- 
resented by short-dashed lines in Fig. 2. The bremsstrahlung 
from these RD and SD clouds (Fig. 2, top-left) fits the data 
in the optical-UV range. The modelling of line spectra in the 
previous section revealed that the dominant component from 
the nuclear region corresponds to RD clouds with V s =100 
kms -1 , no=300 cm~ 3 , and log F h = 10 (short-dash-long- 
dashed lines). 

A small contribution from high velocity clouds with 
V s = 1000 kms _1 and a preshock density no = 1000 cm~ 3 can 
explain the emission in the soft X-ray range (dotted lines). 
This model however does not contribute to the line ratios 
because the line fluxes are low. Moreover, it was found by 
modelling the mid-IR data (Contini et al 2004) from the 
sample of Seyfert galaxies observed by Alonso-Herrero et al. 
(2003) that a high density, low velocity single-cloud model 
(n =1000 cm~ 3 , V s =100 kms -1 ) contributes to the fit of 
most of the galaxies, among them also NGC 4579. This 
model which is also shown in Fig. 2 (long-dashed lines) cor- 
responds to very low line ratios to H/3 and is not included 
in the selected grid (Table 1). 

Recall that the modelling of the spectra have to account 
for covering factors, for the distance of the emitting clouds 
from the AC and for the distance of the galaxy from Earth, 
because the models are calculated at the nebula, while the 
spectra are observed at Earth. The models in Fig. 2 arc 
scaled roughly by taking into account the weights and the 
above mentioned factors. In fact, same types of clouds are lo- 
cated in different regions. Notice that the model correspond- 
ing to the HII regions does not contribute to the continuum 
because its relative weight is low. 

3.2 Infrared bump 

Dust and gas are coupled across the shock front and down- 
stream, so the grains are collisionally heated by the gas to 
temperatures which depend on the shock velocity. Reradi- 
ation from dust in the IR is well explained in Fig. 2 (top- 
right) by the same models which fit the bremsstrahlung. 
For all models a dust-to-gas ratio of 10 -14 by number (cor- 
responding to 4 10~ 4 by mass adopting silicate grains) leads 
to a good fit of the infrared emission bump. In composite 
models corresponding to a relatively low flux from the AC 
or to a low ionization parameter, U, dust grains are mainly 
heated collisionally by the gas throughout the shock front 



and downstream. Therefore, the frequencies corresponding 
to the peaks of the infrared emission depend on the shock 
velocity (see Contini et al. 2004). Thus, it is difficult to dis- 
tinguish between starburst activity and/or an active nucleus 
from the analysis of the infrared bump. 



3.3 Old star contribution 

Fig. 2 (bottom-left panel) shows that a peculiar character- 
istic of NGC 4579 consists in the very high contribution to 
the continuum in the optical range of an old star popula- 
tion background, which is modelled by a black-body corre- 
sponding to T=5000 K (dot-dashed line) . As was explained 
by Contini et al. (2004) Alonso-Herrero et al. data (open 
squares) are uncontaminated from the star contribution. 



3.4 Radio emissions 

The data in the radio range, relative to different apertures, 
are used to determine the different slopes of the contin- 
uum which correspond to the different mechanisms. Fig. 2 
(bottom-right panel) shows a detailed presentation of the 
data and modelling. 

Bremsstrahlung emission from high density gas down- 
stream can be self absorbed at long wavelengths, as was 
found in the Circinus galaxy, NGC 7130, etc. This occurs 
to the downstream emission from clouds with a preshock 
density no ^ 300 cm~ 3 and no = 1000 cm~ 3 . For example, 
model M3 with V s =500 kms~ has a density downstream of 
10 4 cm -3 , a temperature of 10 4 K after recombination in a 
region of about 10 19 cm, leading to optically thick gas (r > 
1) for v sC 10 10 Hz. 

On the other hand, in RD clouds, the bremsstrahlung 
from gas downstream of the low velocity shock explains the 
data from Hummel (1980) and Crane (1977) which refer to 
a region extended to intermediate distances from the AC. 
This is in agreement with Stauffer (1982) who claims that 
in his sample, nuclear radio sources and optical activity are 
well correlated. 

The data relative to the most extended region of the 
galaxy show relatively high fluxes and a different trend. The 
data are well explained by the synchrotron power-law radi- 
ation with a = -0.75 (Bell 1978) (thin solid line) created at 
the shock front by the Fermi mechanism. In fact, the values 
corresponding to large apertures may include anything from 
SN and SNR in the galaxy as a whole. 

The data within 0.1 pc from the AC are investigated 
by Anderson ct al. Although the slope defined by the data 
resembles bremsstrahlung emission with some absorption at 
longer wavelengths, Ulvestad & Ho (2001) claim that the 
free-free interpretation is untenable because the source is 
much more compact (< 0.05 pc) than the dimension of nu- 
clear tori. Advection-dominated accretion flows provide a 
good fit to the spectral slopes but the luminosities are too 
high relative to the black hole masses calculated from the 
diffuse motions. In this paper we adopt the modelling by 
Anderson et al. who explain the observed spectral slope by 
jets (see Sect. 1) 
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Figure 3. The SED obtained by summing up the contribution 
of different models in the different frequency ranges (solid lines). 
The dashed line shows that bremsstrahlung from high density 
clouds is self-absorbed. 

4 DISCUSSION AND CONCLUDING 
REMARKS 

In the previous sections we have modelled the line and con- 
tinuum spectra of NGC 4579. 

Early interpretations of the line spectra (e.g. Stauffer 
1982) based on the high electron temperature derived from 
the low [OIII] 5007/4363 ratio and on the similarity with the 
spectra of SNR, already suggested the presence of shocks. 
However, LINER-like spectra are also shown by galaxies in 
which some filaments, arcs and extended features reveal vig- 
orous recent bursts of star formation (Gonzalez-Delgado & 
Perez 1996, Contini 1997). The weak broad Ha emission ob- 
served in NGC 4579 and UV radiation similar to those of 
Seyfert galaxies suggest an AGN nature of the spectra but 
with stronger low-ionization lines (Barth et al. 2001). 

The results of the present modelling show that the 
multi-cloud models which best explain the line ratios ac- 
count for all types of photoionizing sources, as well as for 
the shocks. The relative importance of the different ionizing 
mechanisms is determined by the relative weight adopted to 
sum up the different single-cloud models. 

The nuclear emission is dominated by the AGN with a 
relative low flux from the AC (log Fh=10) compared with 
that found in the NLR of Seyfert galaxies (log F h ~ 11-13). 
On top of this there is the relative small contribution of a 
young starburst in agreement with Maoz et al (1995) who 
found two emission components in the nucleus of NGC 4579. 
Shock dominated clouds with velocities of 100, 300, and 500 
kms -1 , as well as an older starburst with an age of 4.5 
Myr, probably located in the external nuclear region were 
also revealed by modelling. The good agreement of model 
predictions with the observations is evident in Fig. 4. Fig. 
4 shows a 5" x 5" section of a recent image of the nuclear 




Figure 4. 5" X 5" section of a recent image of the nuclear region 
of NGC 4579, obtained with the Hubble Space Telescope's Ad- 
vanced Camera for Surveys in the F330W band, corresponding to 
~ 3300 A 

region of NGC 4579, obtained with the Hubble Space Tele- 
scope's Advanced Camera for Surveys in the F330W band, 
corresponding to ~ 3300 A (D. Maoz et al., in prepara- 
tion). The extended source 0.6" from the nucleus, previously 
noted by Maoz et al. (1995), is now resolved into a short arc 
which forms the beginning of the ~ 1" spiral arm that can 
be traced for over one revolution around the nucleus. This 
region likely contributes young stellar features to ground- 
based spectra where it is included in the aperture. 

In the extended A and B regions (~ 1 kpc) , RD clouds 
reached by a flux weaker by a factor of 10 than that in the 
nuclear region, due to dilution with distance, contribute to 
most of the line fluxes, while HII regions corresponding to 
T» = 5 10 4 K are also revealed. The contribution of SD clouds 
with V s =100 kms _1 to the [Oil] 3727 and [OIII] 5007 lines 
is significant, confirming that the velocity gradient decreases 
towards the outskirts of the galaxy. Such low velocities are 
generally present in LINERs and LLAGNs (Contini, Viegas, 
& Campos 2003). 

The regions observed by Stauffer (1982) and Keel (1983) 
within intermediate apertures are explained by the same 
AGN corresponding to a low flux, by shocks with a rather 
high V s (500 kms -1 ) and by the HII regions which dom- 
inate the spectra. The high velocity SD clouds created by 
collisional events in SN winds may be related with the strong 
flux from the old stellar population in the optical - near-IR 
range revealed by the continuum SED of NGC 4579. 

Indeed, the calculated V s throughout the galaxy are 
in the range of those observed in the NLR of AGN. Stauf- 
fer (1982) claims that the observed FWHM are 280 kms" 1 , 
while Barth et al find a FWHM of about 520-570 km s -1 for 
the narrow Ha, and for [Nil] and [SB] lines. Pogge et al 
(2000) claim that the bulk of Ha and [OIII] emission comes 
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from a clumpy complex with an elliptical shape of major 
axis of 2" (the emission region is between tens and hun- 
dreds of parsecs). So many different conditions coexist and 
spectra from clouds with different V s can indeed be present. 
The same conclusions were found modelling the Seyfert 2 
galaxies Circinus (Contini, Prieto, & Viegas 1998), NGC 
7130 (Contini et al. 2002), and LINERs (Contini 1997). In 
Seyfert 2 galaxies the AGN dominates in the nuclear region, 
while in LINERs the shock prevails. 

The complex nature of LLAGN as an AGN + starbursts 
is discussed very recently by Cid Fernandez et al. (2004) and 
Gonzalez-Delgado et al. (2004) in the light of new observa- 
tions. They found that few LLAGN have a detectable young 
(< 10 7 yr) starburst component, indicating that massive 
stars do not contribute significantly to the optical contin- 
uum. No features of W-R stars are present. Intermediate 
age populations are very common in LLAGN with relatively 
low [OI]/Hq emission (^ 0.25), but rare for stronger [OI]. 
Moreover, the dominant source of ionization in strong- [OI] 
LLAGN is likely an AGN which is consistent with the de- 
tection of broad Balmer lines in emission in a few cases and 
a larger contribution of older stars in their stellar popula- 
tion. Therefore, it seems that the strong [OI] ([OI]/Ha^ 
0.25) objects are true LLAGN, with stellar processes being 
insignificant. 

Actually, NGC 4579 belongs to the strong [OI] type. 
The models presented in Table 1 which were selected by fit- 
ting the NGC 4579 spectra show that only models Ml and 
M10 correspond to a very low [OI]/Ha. M10 is a matter 
bound model, which does not contribute to NGC 4579 spec- 
tra, and Ml is a SD model corresponding to low V s and no. 
On the other hand, the young starburst corresponding to 
radiation bound model M9, is characterised by [OI]/Ha = 
0.22 (adopting Ha/H/3= 3)), while the power-law domi- 
nated model (M4 in Table 1) and the high velocity shock 
(M3) show high [OI]/Ha line ratios (0.66 and 1.2, respec- 
tively). This confirms that in low [OI]/Ha LLAGN a young 
starburst dominates, as found by Gonzalez-Delgado et al., 
while, starbursts of older ages (t^ 3.3 Myr) show very high 
[OI]/Ha as for AGN. This does not exclude that in com- 
plex high [OI]/Ha objects, a young starburst can be present, 
however, by a low weight, as found for NGC 4579. Also SD 
models show high [OI]/Ha line ratio. These models gener- 
ally explain the emission in supernova remnant, i.e they are 
related with an old stellar population, in agreement with Cid 
Fernandez et al who claim that in the strong- [OI] LLAGN, 
stars cannot play an important role in the gas ionization. 
Recall that the presence of old stars is revealed by the SED 
of NGC 4579. 

We have particularly focused on the nature of the emis- 
sion in the radio range. We have adopted the Ulvestad & Ho 
explaination that synchrotron emission from the base of the 
jet fits the radio data observed within 0.1 pc from the active 
centre. Modelling the continuum SED it is found that radio 
emission within intermediate distances (10-200 pc) from the 
centre is explained by the bremsstrahlung from gas down- 
stream of low velocity shocks (V s =100 kms -1 ) reached by a 
rather low radiation intensity from the AC. The shocks are 
created by the collision of clouds with the jet (Blandford & 
Konigl 1979). Within larger distances from the AC (> 100 
pc) the dominant radio emission is synchrotron radiation 
with a spectral index of -0.75 created by the Fermi mech- 



anism at the shock front, as is often observed from Seyfert 
galaxies (Contini & Viegas 2000). So all types of emission 
observed within different radius from the AC can be rec- 
onciled consistently with the presence of the jet outflowing 
from the accretion disk, if shocks are created by collision of 
the jet with ISM and/or NLR matter. Consistently with the 
presence of the jet, the young starburst can be created in the 
nuclear region close to the base of the jet, by the interaction 
of the high velocity gas with clouds, while the old stars are 
located in the circumnuclear region. The stars were created 
by triggering in collision of shocks with dense matter, so 
the shocks are loosing velocity towards the outskirts of the 
galaxy. 

In conclusion, we have explained both the line and con- 
tinuum spectra in the different regions of the LLAGN NGC 
4579 by composite models accounting for the AGN and for 
the shocks which are created by collision of the jets with ISM 
clouds within different distances from the centre, as well as 
for starbursts of different ages, and for the HII regions. 
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